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Abstract 

The Q-state Potts model in two dimensions in the presence of 
external magnetic fields is studied. For general Q > 3 special choices 
of these magnetic fields produce effective models with smaller Z(Q') 
symmetry (Q' < Q). The phase diagram of these models and their 
critical behaviour are explored by conventional finite-size scaling and 
conformal invariance. The possibility of multicritical behavior, for 
finite values of the symmetry breaking fields, in the cases where Q > 4 
is also analysed. Our results indicate that for effective models with 
Z(Q') symmetry (Q' < 4) the multicritical point occurs at zero field. 
This last result is also corroborated by Monte Carlo simulations. 



1 Introduction 

The ferromagnetic Q-state Potts model in two dimensions is among the most 
studied models of statistical mechanics (see U for a review). In the absence 
of external fields the model has a global Z(Q) invariance [[| which, for low 
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temperatures, is spontaneously broken giving arise to phase trasitions of 
second order for Q < 4 and first order for Q > 4 . The critical fluctuations 
for Q = 2, 3 and 4 are governed by conformal field theories with central 
charges c = |, | and 1, respectively, and the whole operator content of these 
models with several boundary conditions is known [[|. 

In this paper we study the critical behaviour of these models on the 
square lattice, in the presence of symmetry breaking magnetic fields. In 
general these fields will break completely the Z(Q) symmetry, but for some 
special choices of these magnetic fields the resulting effective model will have 
a residual symmetry Z(Q'), with Q' < Q, and a domain wall structure at 
low temperatures, similar to these in the Q'-state Potts model. 

On general grounds |5[], we do expect that in order-disorder phase transi- 
tions, of discrete symmetry models, the critical behavior is dictated mainly 
by the number of ground-states (zero temperature configurations) and the 
relative surface energy of the infinite domain walls connecting these ground 
states. This reasoning induces us to expect, for arbitrary values of the sym- 
metry breaking fields, producing a Z(Q') model, an effectice model in the 
same universality class as the Q'-state Potts model. However this analysis 
is not valid in general since some critical and multicritical models, like the 
Ising model and the tricritical Ising model ||, although having the same 
number of ground-states and domain wall structure at low temperatures, 
exhibit distinct critical properties. 

An earlier mean field analysis on these models @ indicate that for Q < 4 
the Z(Q') (Q' < Q) model, produced by the breaking fields, are in the same 
universality class of the Q'-state Potts model, for nonzero values of the fields, 
while for Q > 4, with Q' < 4, there is a multicritical point for finite values 
of the critical field, where the phase transition changes from first to second 
order. 

Our study will be done numerically by using standard finite-size scaling 
[H] to obtain the phase diagram of the models, and the machinery arising 
from conformal invariance || |10| to distinguish the several possible critical 



behaviours. In the location of multicritical points for Q > 4 we also perform 
Monte Carlo some simulations by calculating the fourth-order cumulant of 
the magnetization. 
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2 The model 



Defining at each lattice site f — of a square lattice an integer variable 
nff = 0,1, ... ,Q — 1, the Hamiltonian of the Q-state Potts model with nh 
(nh — 0, 1, . . . , Q — 1) symmetry breaking fields {h m } (m = 0, 1, . . . , rih — 1) 
is given by 

H Q (e, {h m }) = -e ^ <W,rv - £ £ h m<W,m- (1) 

<f,f'> m=0 t 

where e > is the ferromagnetic coupling and the first sum runs over nearest- 
neighbour sites. In the absence of external fields the model has a Z(Q) sym- 
metry, since the configuration {n?} and {n? + 1, mod.Q} (I = 1, 2, . . . , Q) has 
the same energy. The fields {h m }, depending on their relative values, break 
this symmetry totally or partially. The interesting cases where the symme- 
try is partially broken are those where h m =h > (m = 0, 1, . . . , — 1) and 
the remaining symmetry is Z^n^) ® Z(Q — n/J. This symmetry corresponds 
to Zirih) rotations among the variables pointing in the field directions and 
Z(Q — rih) rotations among the other variables. At zero temperature we have 
rih ground states and we do expect that the Z{rih) symmetry is spontaneously 
broken. 

Rather than working with the above Euclidean version of the model it is 
convenient to consider its quantum Hamiltonian version in order to simplify 
our numerical analysis. The row-to-row transfer matrix as well the associ- 



ated r-continum quantum Hamiltonian [ I I ] can be derived by a standard 



procedure (see |12| for example). The associated one- dimensional quantum 



Hamiltonian in a L-site chain is given by 

L q-l , n h -l . . 

#=-££((^ + i) a + A£f+ E h m (Sie-^ m ) a ), (2) 

1=1 a=0 ^ m=0 ' 

where A plays the role of temperature and the magnetic fields {h m } (m = 
0,1, . . . ,n h — 1) are related with the fields {h m } in ([!]). In (|2|) 5*; and Ri are 
Lft ® Lft matrices satisfying the Z(Q) algebra 

[Ri, St] = (@ — ^)Sk,iSkRi Ri = Sf = 1. 

where 9 = exp(^), and in the basis where Si is diagonal they are given by 

Si=l®l®---I®5 , ®l---®1, 
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i2j=l®l<g>...l<g>.R®l---<g>l 

where the matrices S and R are in the I th position in the product and are 
given by 

S = E&0 1 0% ><i\, R = EgTo 1 \i >< [i + 1]q|, 

and the symbol [x + u]q means the addition {x + y), modulo Q. 

In the absence of magnetic fields (h m = 0; m = 0, 1, . . . , n h — 1) the Z(Q) 
symmetry of (|T|) is reflected in (EI) by its commutation with the Z(Q)-charge 
operator 

P = UtiRi, P Q = ^- (3) 

The Hilbert space associated with (EJ) can therefore be separated into disjoint 
sectors labelled by the eigenvalues exp(^=p), (q — 0, 1, . . . , Q — 1) of (^). The 
interesting cases, which we will concentrate on in this paper, are obtained by 



choosing in (g) equal values for the magnetic fields [13[ 

h x = h 2 = . . . = K h = h > 0. (4) 

In this case the symmetry Z{n h ) (g) Z(Q — n h ) of (|I|) is reflected by the 
simultaneous commutation of (fj) with the "parity" operators 

V = Ui=iV h V n » = l, (5) 

w = nf=i w h w Q - nh = i, (6) 

with 

Vj = 1®1®---1®V®1---®1 

W", = l®l®...l<g>W<g>l---<g>l 

and V and W, located at the Z t/l position in the product, are Q x Q matrices 
given by 

n h -l Q-l 

V= l*>< [^ + 1 ]"J+ E |i><*l 

i=0 i=n-h 



nh— 1 Q-l 

^ U" N» i" I _l_ \ U' ^ ^ U _L 1 1 „ 



VP = ^ |i >< i\ + \i >< [* + 1] 

i=0 i =n h 
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The Hilbert space associated to (0) is now separated into Uh(Q — Uh) disjoint 
sectors labelled by the eigenvalues exp(^^) and exp(^p^-) (v = 0,1, ... ,rih— 

1, w = 0, 1, . . . , Q — rih — 1 ) of the operators V and W , respectively. 

In the numerical diagonalization of (0), with periodic boundaries, all the 
above symmetries, together with the translational invariance, enables us to 
handle large lattices with modest computer time and memory. We use the 
Lanczos method to diagonalize (FJ) up to L = 10, 11 and 13, for Q = 5, 4 and 
3, respectively. 

3 Results 

We considered in our study only the interesting cases where Q > nh > 2 
and hi = h 2 = . . . = h nh — h > 0, since in these cases always a remaining 
symmetry Z(Q') (Q' = Uh > 2) still remains for h ^ 0. 

When h = the model is self-dual with a phase transition at A c (0) = 1, 
which has a second order or first order nature depending if Q < 4 or Q > 4, 
repectively. In the limit h — > oo the eigenvectors of Si in (0) with eigenvalues 
9 l , I = rih, rih + 1, . . . , Q — 1, are forbidden and we have an effective n^-state 
Potts model at zero field. Analysing the effect of (0) in the remaining Tri- 
dimensional Hilbert space it is not difficult to see that the phase transition 
happens at 

X c (h -00) = -^ (7) 
n h 

Between those two extremum values of h we estimate the phase transition 
curve X c (h) by using standard finite-size scaling. The curve is evaluated by 
extrapolations to the bulk limit (L — ► 00) of sequences X c (h, L) obtained by 
solving p| 

r L (A c )L = r L+1 (A c )(L + l), L = 2,3..., (8) 

where r^(A c ) is the mass gap of the Hamiltonian (^|) with L sites. 

Once the transition curve is estimated, in the region of continuous phase 
transitions we expect the model is conformally invariant. This symmetry 
allows us to infer the critical properties from the finite-size corrections to the 
eigeinspectrum at A c ||. The conformal anomaly c can be calculated from 
the large L behaviour of the ground-state energy Eq(L). For periodic chains 
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Eq(L) behaves as 

— = — + ), (9) 

where e M is the ground-state energy, per site, in the bulk limit and v s is the 
sound velocity. The scaling dimensions of operators governing the critical 
fluctuations (related to critical exponents) are evaluated from the finite-!/ 
corrections of the excited states. For each primary operator, with dimension 
x,/, and spin s^, in the operator algebra of the system, there exists an infinite 
tower of eigenstates of the quantum Hamiltonian, whose energy E^ m , and 
momentum P^ m /, in a periodic chain are given by 



< m ,(L) = E + ^{ X<t> + m + m>)+ o{L^) 

P m,m' = ( S 4> + 171 ~ m ') 



(10) 



where m, m' — 0, 1, . . . . 

We present our results separately for the cases Q < 4 and Q > 4 in the 
next sections. 



3.1 Models with Q < 4 

There exist three interesting cases, namely, the 3-state Potts model with two 
fields (Q = 3, rih = 2) and the 4-state Potts model with three and two fields 
(Q = 4,n h = 3,2). 

The critical curves were obtained by solving <^). As an example, in Fig.l 
we show the extrapolated curve for the case of Q = 3 and n h = 2. We also 
show in this figure the curve obtained by solving (^) for L = 5. We clearly 
see an agreement with the limiting values A c (0) = 1 and X c (h — > oo) 
predict by (0). Similar curves are obtained in the other cases. 

The conformal anomaly and anomalous dimensions are obtained using 
relations (JS , ITI ) , for several values of h. In table 1,2 and 3 we show for some 
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values of h the extrapolated results obtained for the cases (Q — 3,71^ = 2) 
(Q = 4,rih = 3) and (Q = 4, = 2), respectively. Our personal estimative 
of errors are in the last digit. In these tables we also present our conjectured 
values. The dimensions x n (k,v) and x n (k,v,w) appearing in these tables 
are obtained by using in ( PH|) the nth (n = 1,2,3,...) eigenenergy in the 
sector with momentum ^ (k = 0, 1 . . .) and eigenvalues exp(^p^-) and 
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exp( !?™ )(v = 0, 1, . . . , rih — 1, w = 0, 1, . . . , Q — rih — 1) of the operators V 

and W denned in (|5|) and (||), respectively. 

We see in table 1 that for all values of h the conformal anomaly is c = |, 
indicating that the model share the same universality class as the Z(2) Ising 
model. The dimensions 1 and | correspond to the dimensions of the energy 
and magnetic operator in the Ising model and the dimension 2=1+1 and 
| = 1 + | are the next dimensions in the tower of these operators (see Eq. 
©)■ 

In table 2 we clearly see that the conformal anomaly, for all values of 
h is c = |. There exists two modular invariant universality classes of con- 
formal theories with c = | |L4|]. One of these can be represent by the re- 
stricted solid-on-solid (RSOS) model [|15|, [16| and the other by the 3-state 
Potts model with no magnetic fields. These two models, although having 
the same conformal anomaly c = |, have distinct operator content. The di- 
mension x = 0, |, y , |, | in table 2 and the degeneracy of the sectors where 
the operator V in (||) has eigenvalues exp(y E ) and exp(y L ) indicate that the 
model belongs to the same universality class as the 3-state Potts model. 

In table 3, as in table 1, the conformal anomaly is c = 1/2 and the 
dimensions are those of the Ising model indicating that both models are in 
the same universality class. 

Beyond the values of h presented in tables 1-3 we also performed a care- 
full analysis for small values of (h ~ 0.01) finding similar results as those 
presented in those tables. This indicate that for Q < 4, where the model has 
a second-order phase transition in the absence of external fields, the intro- 
duction of rih (Q > Uh > 2) fields, of arbitrary strength, brings the model 
into the universality class of a Potts model with rih states. 



3.2 Models with Q > 4 

In this case, while in the absence of external fields the models exhibit a first- 
order phase transition the introduction of magnetic fields (4 > > 2) 
of infinite and equal strength (h — > oo) render them to an effective n^-state 
Potts model, which has a second order phase transition. This brings the 
interesting possibility of a multicritical behaviour for a finite value of h, 
when the transition curve changes from second to first order as we decreases 
h from the infinite value. In fact this is the mean field prediction 0. 
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Since the Hilbert space grows exponentially with Q, the simplest case 
where the above critical point may occur is the 5-state Potts model in the 
presence of = 2 magnetic fields. In table 4 we present, for some values of 
h, our results for the finite-size sequences of the conformal anomaly of the 
model. These sequences are obtained from ([U]) and @ 

12[(L+1)E (L)-LE (L + 1)] 
L ' L+l [(L + iy-L 2 }[E 2 (L + l)-E 1 (L + l)} { ' 

where Eq(L) is the ground-state energy for the chain of length L and Ez(L), 
Ei(L) are the lowest eigenergies with momentum and respectively, in the 
sector where the operators V and W, defined in (|5]) and (^[) has eigenvalues 
(-1) and (1). We see in table 4 that for h > 0.05 we still have an Ising-like 
behaviour with c = ~. For h < 0.05, and for the lattice sizes we were able to 



handle, it is not possible to obtain reliable results using (|iT|). 

An heuristic method, which was proved to be effective in obtaining mul- 
ticritical points in earlier works [1% is to simultaneously solve Eq. (||) for 



three different lattice sizes 

r L (A c )L = r L+1 (A c )(L + 1) = r L+2 (A c )(L + 2). (12) 

We tried to solve these equations for 0.5 > A > (L = 5) and we found no 
consistent solutions, which indicates the absence of a tricritical point for a 
finite value of h. 



Another method, also used to locate multicritical points [17], is obtained 
from the simultaneous crossing of two different gaps on a given pair of lattices 
(instead of three lattices as in (|T2"D). Trying several different gaps we also 
did not find, within this method, a multicritical point for h ^ 0. 

Since these methods are heuristic and the lattice sizes we are considering 
may not be enough to obtain the bulk limit (L — > oo) in the region h ~ 
we decide to supplement our results by Monte Carlo simulations. These 
simulations will enable us to distinguish the order of the phase transition as 
we change the magnetic field strength. We simulate the systems by the heath- 
bath algorithm and analyse the fourth-order cumulant of the magnetization 
as a function of the magnetic field. The simulations were done directly in the 
classical version of the model (|I|). Since the evidence of a multicritical point, 
for non-zero values of h, would be large for higher values of Q, we choose 
Q = 7 and rih = 2 (ho=hi = h) for extensive calculations. 
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The fourth-order cumulant of the magnetization is defined by 

where the averages are done on an L x L lattice. The magnetization in 
for a given configuration {n,?} of classical variables, is defined by 



Following Binder |T8[ , the cumulant (|T3| ) will be zero for T > T c and Vl — \ 
for T < T c . At the transition temperature T c (|T3|) will be zero for continous 
phase transitions and negative for first-order phase transitions. 

In Fig. 2 we show the values of Ul for lattice sizes 50 x 50, 80 x 80 and 
170 x 170. In the simulations we choose e = 1, h=0.01 and each point was 
obtained by averaging 5.10 iterations, after thermalization. We see in this 
figure that while for L = 50 the phase transition appears to be first order, as 
L grows the numerical results indicates the phase transition to be continuous. 
The result for the smaller lattice L = 50 is clearly due to the finite size of the 
lattice. By repeating these simulations for even smaller values of h we should 
expect that these finite-size effects will be apparent for even larger lattices, 
and our simulations are in favour of a multicritical point only at h=0. 



4 Conclusion 

We have calculated the phase transition diagram and critical properties of the 
Q-state Potts model in the presence of (Q > > 1) external magnetic 
fields of equal strengh h > 0. In the case where Q > rih > 2 the original 
symmetry, at h = 0, breaks into a Z{rih) <8> Z(Q — rih) symmetry. The Z(rih) 
part of the above symmetry relates the configurations of the distinct 
ground-state configurations at zero temperature, and by standard arguments, 
should be spontaneously broken at low temperatures. 

Our results, based on conformal invariance and supplemented by Monte 
Carlo simulation indicate that, for arbitrary values of h, the order- disorder 
phase transition associated with the global Z(nh) symmetry is in the same 
universality class of the n^-state Potts model. Morever, for Q > 4, contrary 
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to the mean field prediction, we do not see any evidence of a multicritical 
point for non-zero values of h. 
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Figure Captions 



Fig. 1 - Estimates for the critical curve of the 3-state Hamiltonian (|^) 
with rih = 2 fields. The curve in the largest scale, for < h < 5, interpolates 
the points obtained by extrapolating the solutions of (^) for L = 2 — 13 
(circles). The inserted curve for < h < 50 was obtained by solving (||) for 
L = 5. 



Fig. 2 - Fourth-order cumulant of the magnetization (|13| ) as a function of 
(3 for the 7-state Potts model in presence of rih = 2 external fields ho=hi=0.01 
(seeEq. (|)). The lattice sizes are L = 50x50, L = 80x80 and L = 170 x 170. 
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Table Caption 



Table 1 - Extrapolated and conjectured results for the conformal anomaly 
c and anomalous dimensions x n (k,v) of the 3-states Potts chain (|2|) with 
n h = 2 magnetic fields (h — h\ — h) (see the text). The conjectured values, 
in parenthesis, are the corresponding ones in the critical Ising model. 



Table 2 - Extrapolated and conjectured results for the conformal anomaly 
c and anomalous dimensions x n (k,v) of the 4-states Potts chain (|2|) with 
n h = 3 magnetic fields (h — hi — h 2 = h) (see the text). The conjectured 
values, in parenthesis, are the corresponding ones appearing in the 3-state 
Potts model. 



Table 3 - Extrapolated and conjectured results for the conformal anomaly 
c and anomalous dimensions x n (k,v,w) of the 4-states Potts chain (0) with 
n h = 2 magnetic fields (h — hi = h) (see the text). The conjectured values, 
in parenthesis, are the corresponding ones in the critical Ising model. 



Table 4 - Finite-size sequences cl,l+i, defined in (|j) for the 5-state Potts 
model Hamiltonian (§) with rih = 2 (ho = hi = h) external magnetic fields. 
The last line gives the extrapolated results. 
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Table 1 



h 


c 


x 2 (0,0) 


xi(l,0) 


xi(0,l) 


xi(l,l) 


0.5 


0.50087 


1.00073 


2.001 


0.1251 


1.1251 




(0.5) 


(1) 


(2) 


(0.125) 


(1.125) 


1.0 


0.50000 


0.99996 


1.9995 


0.12497 


1.12497 




(0.5) 


(1) 


(2) 


(0.125) 


(1.125) 


1.5 


0.50003 


1.00000 


1.99997 


0.1249 


1.1249 




(0.5) 


(1) 


(2) 


(0.125) 


(1.125) 


2.0 


0.50008 


1.00000 


2.00000 


0.12500 


1.12500 




(0.5) 


(1) 


(2) 


(0.125) 


(1.125) 



Table 2 



h 


c 


^(0,0) 


^(0,0) 


Z2(0,l) 


2:3(0,1) 


2:i(l,l) 


2:2(1,1) 


0.5 


0.800 
(0.8) 


0.8001 
(0.8) 


2.79 
(2.8) 


0.1333 
(0.133...) 


1.332 
(1.33...) 


1.1333 
(1.133...) 


2.332 
(2.33...) 


1.0 


0.7991 
(0.8) 


0.7997 
(0.8) 


2.8 
(2.8) 


0.1329 
(0.133...) 


1.330 
(1.33...) 


1.1329 
(1.133...) 


2.331 
(2.33...) 


1.5 


0.79 
(0.8) 


0.799 
(0.8) 


2.83 
(2.8) 


0.1334 
(0.133...) 


1.330 
(1.33...) 


1.1334 
(1.133...) 


2.3319 
(2.33...) 


2.0 


0.7991 
(0.8) 


0.7999 
(0.8) 


2.84 
(2) 


0.13345 
(0.133...) 


1.330 
(1.33...) 


1.13345 
(1.133...) 


2.3319 
(2.33...) 
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Table 3 



h 


c 


x 2 (0,0,0) 


m(l,0,0) 


m(0,l,0) 


xi(l,l,0) 


0.5 


0.5002 
(0.5) 


0.997 
(1) 


1.944 

(2) 


0.1248 
(0.125) 


1.1248 
(1.125) 


1.0 


0.505 
(0.5) 


1.0000 

(1) 


1.999 
(2) 


0.12509 
(0.125) 


1.12509 
(1.125) 


1.5 


0.5000 
(0.5) 


0.9999 
(1) 


1.9996 
(2) 


0.1249 
(0.125) 


1.1249 
(1.125) 


2.0 


0.5000 
(0.5) 


0.9999 
(1) 


1.9997 
(2) 


0.12499 
(0.125) 


1.12499 
(1.125) 



Table 4 





0.05 


0.1 


0.5 


1.0 


4 


1.494578 


1.293239 


0.822347 


0.747102 


5 


1.150201 


0.947010 


0.665147 


0.632826 


6 


0.957555 


0.784083 


0.603697 


0.583638 


7 


0.834717 


0.696948 


0.572906 


0.558058 


8 


0.753559 


0.646537 


0.554626 


0.542796 


9 


0.698775 


0.615352 


0.542598 


0.532888 


L 10 


0.660931 


0.594802 


0.534149 


0.526071 


oo 


0.5693 


0.5448 


0.4898 


0.500 
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